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During  the  course  of  this  NRL/DARPA  program,  we  developed  a  novel,  kinetic  Monte 
Carlo  (kMC)  simulation  code  capable  of  modeling  diamond  growth  on  the  atomic  scale.  We  also 
developed  a  two-dimensional  simulation  code  capable  of  predicting  the  evolution  of  faceted 
polycrystalline  diamond  films  as  a  function  of  diamond  growth  conditions. 

We  developed  a  kinetic  Monte  Carlo  technique  for  simulating  thin  film  growth  by  chemical 
vapor  deposition  which  is  inherently  atomistic  and  three  dimensional,  and  incorporates  known 
reaction  kinetics.  The  kinetics  of  the  relevant  surface  processes  and  the  effects  of  surface  atomic 
structure  and  morphology  are  included  within  this  model.  Processes  which  involve  cooperation 
between  adjacent  surface  atoms,  such  as  reconstruction  and  growth  of  { 100}  faces  of  diamond,  are 
handled  explicitly  and  atomistically.  The  temporal  evolution  of  the  system  is  simulated  using  an  N- 
Fold  Way  Monte  Carlo  scheme  which  employs  a  variable  time  increment,  allowing  the 
consideration  of  processes  that  occur  on  widely  different  time  scales.  The  model  is  efficient  enough 
to  handle  the  growth  of  hundreds  of  atomic  planes  on  substrates  containing  thousands  of  surface 
atoms,  using  only  a  desktop  workstation.  Because  the  model  is  inherently  atomistic,  it  is  useful  not 
only  for  studying  the  evolution  of  film  morphology  and  defect  incorporation,  but  also  as  a  tool  for 
examining  and  verifying  proposed  atomic  growth  mechanisms.  The  growth  of  {111  {-oriented 
diamond  films  was  simulated  for  several  substrate  temperatures  ranging  from  800  to  1500  K. 
Under  typical  CVD  conditions,  the  growth  rates  vary  from  about  0.1  to  1.2  pm/hr  between  800 
and  1500  K,  and  the  activation  energy  for  growth  on  the  {111}  surface  between  800  and  1 100  K 
is  11.2  kcal/mole,  in  agreement  with  experiment  and  other  simulations.  The  concentrations  of 
incorporated  point  defects  are  low  at  substrate  temperatures  below  1300  K,  but  become  high  above 
this  temperature.  While  the  growth  rates  are  controlled  by  the  combined  contributions  to  growth 
from  CH3  and  C2H2,  the  incorporation  of  point  defects  depends  on  competition  between  C2H2  and 
CH3  growth.  Therefore,  in  the  growth  environment  simulated  here,  optimal  {111}  growth 
efficiency,  defined  as  the  ratio  of  growth  rate  to  point  defect  concentration,  can  be  achieved 
between  approximately  1100  and  1200  K. 

The  growth  of  {001 },  {Oil },  and  {111}  diamond  films  was  simulated  in  various  gaseous 
environments  containing  H,  H2,  CH3,  and  C2H2  concentrations  typical  of  hot  filament  CVD 
reactors.  The  effects  of  surface  chemistry  is  incorporated  in  the  model  through  conventional 
reaction  rate  coefficients,  and  surface  atomic  structure  and  morphology  are  included  by 
representing  the  films  atomistically  in  three  dimensions.  The  formation  of  dimer  bonds  between 
surface  atoms  is  handled  explicitly,  allowing  arbitrary  dimer  surface  reconstructions  of  {001 } 
facets.  The  concentrations  of  H,  CH3,  and  C2H2  were  varied  independently  over  six  decades,  and 
the  effects  of  these  variations  on  growth  rate  and  film  morphology  were  ascertained.  Changes  in  H 
concentration  affect  growth  rates  the  most,  since  the  rate  of  formation  of  radical  sites  on  the 
diamond  surface  is  proportional  to  H  concentration.  Rates  of  growth  of  all  facets  are  comparable  at 
each  gas  composition,  except  when  the  concentration  of  C2H2  is  low  (below  approximately  1x10" 


2%).  The  disparity  in  growth  rates  at  low  C2H2  concentrations  is  due  to  the  coordination  that 
adsorbed  C  atoms  must  assume  to  form  a  diamond  cluster  on  a  flat  surface.  The  {001 },  {01 1 }, 
and  {111}  films  require  that  one,  two,  and  three  neighboring  adsorbed  C  atoms,  respectively, 
bond  together  in  order io-initiate  a  monolayer  on  a  flat  face.  However,  the  chemisorption  of  C2H2 
contributes  two  carbon  atoms,  and  thus  this  effect  is  only  significant  at  low  C2H2  concentrations. 
The  growth  of  {111}  films  in  C2H2-deficient  atmospheres  is  controlled  by  the  nucleation  of 
diamond  clusters  on  atomically  smooth  faces.  Thus  {111}  films  grow  by  step  flow  under  these 
conditions,  and  are  atomically  smooth  for  extended  periods  (up  to  several  minutes)  during  growth. 

The  growth  of  diamond  films  and  the  rate  of  incorporation  of  vacancies,  H  atoms,  and  sp2 
defects  were  simulated  for  several  substrate  temperatures  and  feed  gas  compositions.  An  efficient 
atomic-scale  on-lattice  kinetic  Monte  Carlo  method  was  used  to  capture  a  simple  yet  physically 
meaningful  representation  of  the  atomic  processes  involved  in  defect  trapping.  Growth  on 
(100)(2xl)-  and  (lll)-oriented  diamond  surfaces  was  simulated  for  substrate  temperatures  from 
600  °C  to  1200  °C  with  0.4%  inlet  CH4,  and  for  0.4%  to  7%  inlet  CH4  concentrations  with  a 
substrate  temperature  of  800  °C.  The  simulations  predict  that  H  atom  concentrations  increase  with 
substrate  temperature  as  the  relative  contribution  of  CjHj  to  deposition  increases,  simply  because 
the  larger  QHj  molecule  can  cover  defect  sites  more  easily  than  can  CH3.  The  concentrations  of 
incorporated  sp2  bonds  increase  with  the  CH4  content  in  the  feed  gas  because  the  concentration  of 
H,  which  is  required  to  convert  sp2  bonds  to  sp3  bonds  at  the  surface,  decreases  with  increasing 
inlet  CH4.  Vacancy  concentrations  are  low  under  all  simulated  conditions.  The  ratio  of  growth  rate 
to  H  atom  concentration  is  highest  around  800-900  °C,  and  the  growth  rate  to  sp2  ratio  is  maximum 
around  1%  CH4,  suggesting  that  these  conditions  are  ideal  for  diamond  growth  under  the  simulated 
conditions,  as  commonly  observed  experimentally. 

Chemically  vapor  deposited  (lOO)-oriented  diamond  is  known  to  grow  more  slowly  than 
(1 10)-oriented  material  and  to  contain  smooth  surface  features.  The  existing  self-consistent  models 
for  diamond  growth,  on  the  other  hand,  predict  that  (lOO)-oriented  films  grow  fastest  and  are 
atomically  rough.  Surfaces  which  grow  slowly  and  develop  smooth  features,  such  as  (111),  do  so 
because  the  nucleation  of  new  growth  layers  on  these  surfaces  is  much  slower  than  growth  at  step 
edges, but  there  is  no  nucleation  barrier  for  the  formation  of  new  growth  layers  on  (100)  surfaces. 
This  suggests  that  the  experimentally  observed  (100)  growth  behavior  might  be  rationalized  by  the 
introduction  of  growth  mechanisms  that  hinder  layer  nucleation  and  promote  step-flow  growth. 

One  such  mechanism  is  the  preferential  etching  of  monomolecular  C  islands.  The  thermochemistry 
of  the  (100)  dimer  insertion  reaction,  by  which  CH2  species  are  incorporated  at  the  surface, 

indicates  that  the  incorporation  reaction  is  reversible  and  that  an  etching  mechanism  might  indeed 
exist.  We  have  presented  the  results  of  coupled  tight  binding  (TB)  and  density  functional  theory 
(DFT)  quantum  mechanical  calculations  that  suggest  that  etching  of  isolated,  monomolecular  C 
islands  occurs  at  an  appreciable  rate,  while  etching  from  larger  C  islands  (i.e.,  step  edges  and 


terraces)  is  not  favorable.  The  TB/DFT  calculations  and  the  proposed  etching  mechanism  were 
combined  with  previously  determined  hydrocarbon  reaction  rates  on  diamond  surfaces  and  input  to 
a  kinetic  Monte  Carlo  model  of  the  chemical  vapor  deposition  of  diamond  on  (100),  (1 10),  and 
(1 1 1)  surfaces.  The  simulation  results  clearly  demonstrate  that  the  inclusion  of  etching  produces 
much  smoother  (100)  diamond  surfaces  (as  compared  with  calculations  in  which  etching  is 
omitted).  Etching  slows  (100)  growth  considerably  without  substantially  affecting  (110)  and  (111) 
growth  rates.  The  growth  rates  that  are  predicted  when  etching  is  included  in  the  model  are  in  very 
good  agreement  with  a  wide  range  of  experimental  evidence  that  suggest  that  (100)  is  slower  than 
(110)  growth  and  comparable  to  (1 1 1)  growth  under  hot  filament  CVD  reactor  conditions. 
However,  while  etching  is  key  to  our  understanding  of  diamond  growth  in  general,  and  (100) 
growth  in  particular,  several  details  of  (100)  growth  remain  unexplained. 

We  performed  a  series  of  simulations  of  the  growth  of  polycrystalline,  faceted  films  from 
randomly  oriented  nuclei  in  two  spatial  dimensions.  The  simulations  track  the  motion  of  all  comers 
where  facets  from  the  same  grain  and  different  grains  meet.  While  earlier  studies  have  employed 
this  type  of  simulation  method,  this  study  provides  the  most  complete  analysis  of  the  resultant 
microstructures  performed  to  date.  Examination  of  the  microstructures  showed  that  the  grain 
boundaries  are  straight,  but  may  turn  abruptly  when  facet  comers  meet,  and  grains  are  annihilated 
when  two  grain  boundaries  meet  (an  observation  made  earlier  by  van  der  Drift).  Much  of  the 

growth  behavior  may  be  understood  in  terms  of  the  a-dependent  idiomorphs  (the  shape  found  by 
convexifying  a  polar  plot  of  growth  velocity  versus  surface  orientation).  The  surface  morphology 
is  sensitive  to  the  value  of  a,  especially  for  a’s  where  idiomorphs  suffer  a  change  in  topology. 

The  mean  grain  size  and  RMS  surface  roughness  were  found  to  grow  as  4h  (where  h  is  the  film 
thickness)  in  agreement  with  experiment,  theoretical  and  previous  simulation  results.  The  grain  size 
distribution  was  found  to  be  temporally  self-similar  when  the  grain  sizes  were  scaled  by  the  mean 
grain  size.  The  grain  size  distribution  is  better  fit  by  a  gamma  distribution  than  by  a  log-normal 

distribution.  The  crystallographic  orientation  distribution  (i.e.,  texture)  peaks  at  an  a-dependent 
orientation  and  the  peak  sharpens  during  film  growth.  The  orientation  distribution  is  well  described 
by  a  Gaussian,  in  agreement  with  theory.  The  a-dependent  peak  position  may  be  derived  directly 

from  the  orientation  of  the  largest  radius  vector  of  the  appropriate  idiomorph  in  two  and  three 
dimensions.  The  volume  of  material  grown  from  each  facet  type  exhibits  a  cross-over  at  finite 

thickness  to  that  predicted  on  the  basis  of  the  a-dependent  idiomorphs,  in  agreement  with 
experiments. 
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